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Abstract - We study a generalized Thomson problem that appears in several condensed matter 
settings: identical point-charge particles can penetrate inside a homogeneously charged sphere, 
with global electro-neutrality. The emphasis is on scaling laws at large Coulombic couplings, 
and deviations from mean-field behaviour, by a combination of Monte Carlo simulations and 
an analytical treatment within a quasi-localized charge approximation, which provides reliable 
predictions. We also uncover a local overcharging phenomenon driven by ionic correlations alone. 
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Ph ' The venerable Thomson problem of finding the ground 
^ state of an ensemble of electrons confined in an homoge- 
i— heously charged neutralizing sphere, is still unsolved and 
has a long history, see e.g. [IH3] and references therein for 
^.the different generalizations that have been put forward. 
^^Thc model was introduced at the beginning of the 20th 
\Q> century [T], just after the discovery of the electron, but 
before that of the proton, as a classical representation of 
C*"Hhe atom; hence the "plums" representing the electrons, 
"xj- and the introduction of an homogeneous background (the 
l O ''pudding"), to fulfill electro- neutrality. This picture, al- 
though obsolete in the atomic realm, has nevertheless at- 
. . tracted the interest of mathematicians, physicists and bi- 
K^qlogists alike, due to its relevance in particular for ionic 
K> ordering at interfaces [I] , for the behaviour of colloids self- 
assembled at the edge of emulsion droplets (colloidosomes, 
see e.g. [5j), for the study of one-component plasmas and 
their experimental realizations (electrons on a liquid He- 
lium surface [6]) or for understanding viral morphology 
[7]. The Thomson problem reappeared recently in sheep's 
clothing in different contexts, from the screening effects 
in hydrophobic polyelectrolytes [Hj, to the behaviour of 
Coulomb balls (identical particles confined in an harmonic 
trap [TTj). including hydrogels [15] . where the uptake of 
counter-ions by a cross-linked polymer network (the "pud- 
ding") is the key feature leading to the expansion of the 
network by osmotic pressure, hence the capability to ab- 
sorb large quantities of water. At variance with Thomson's 



preoccupation where ground state configurations were un- 
der scrutiny, those articles were concerned with the finite 
temperature behaviour (T ^ 0). At high to moderate tem- 
peratures, mean-field theory provides a trustworthy frame- 
work, and allows to obtain some analytical results [S]. The 
regime of low temperatures (large couplings, a notion to be 
specified below) is more elusive, and it will be our primary 
objective in the present contribution. Several analytical 
predictions will be derived, including scaling laws for two 
important quantities characterizing the screening proper- 
ties. These predictions will be tested against numerical 
simulations, that also give access to detailed microscopic 
information concerning the structure. 

We start by defining the model and introducing a rel- 
evant coupling parameter. We consider a single perme- 
able and spherical globule of radius R g and charge —zqZ g 
(referred to as the background), surrounded by its Z g 
counter-ions of charge zq, where q is the elementary charge 
and z the ionic valency. It should be stressed that the 
counter-ions can penetrate but also leave the homoge- 
neously charged globule, which is an important difference 
with Thomson's original formulation. Without loss of gen- 
erality, the globule is assumed negatively charged (positive 
counter- ions). In such a salt-free system, a confining cell 
is required to avoid evaporation of all counter-ions; its ra- 
dius is denoted R c . A key quantity is the "globule total 
uptake" charge — Z up , which includes the background and 
the counter-ions present within the globule [13]. Hence, 
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Z g — Z up can be viewed as the number of counter-ions out- 
side the globule. Following [TUIITT] . we define the coupling 
(plasma) parameter T as the ratio between the character- 
istic electrostatic energy of a counterion-counterion inter- 
action, z 2 q 2 /(4ire5), and the thermal energy, kT . Here 
e is the solvent dielectric constant, and the distance 8 is 
taken as the ion sphere radius: Z g S 3 = i? 3 . Introducing 
the Bjcrrum length Is = q 2 /(47refcT), which is about 7 A 
in water at room temperature, we get 



*V = -1/3 
4Tre5kT 9 R g ' 



(1) 



For r < 1, mean- field Poisson-Boltzmann theory [10] pro- 
vides an accurate description, while the strong coupling 
regime corresponds to T > 1. It has been shown in [5] 
that within the non- linear mean- field regime [5], counter- 
ion penetration into the globule lead to Z up oc \[Z~ g for 
large enough bare charge Z g . For the purpose of compar- 
ison with numerical calculations in the strongly coupled 
regime, this quantity needs suitable rescaling, and we de- 
fine 



n = z uv 
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which only depends on T [TJ]. Within mean-field, we then 
have II oc r 3 / 2 up to a prefactor of order 1, and in addition, 
the relation n oc exp(—ze(j)/kT) between counter-ions den- 
sity n(r) and the local mean electrostatic potential 4>{r) 
[TO] , allows us to relate the reduced charge II to the char- 
acteristic decay length of n(r) at r 
define the positive quantity 



R g . To this end, we 
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log n(r), 
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that can be viewed as the reduced (inverse) decay length of 
the counter-ion profile in the globule vicinity. Gauss theo- 
rem implies S = II, again within irrelevant prefactors: the 
reduced decay length S" -1 is inversely proportional to the 
reduced total charge, that is itself an increasing function 
of the background charge Z g . 

How are the previous results affected in the strong- 
coupling regime? To make analytical progress when T 3> 
1, we take advantage of the caging of particles that takes 
place under their strong mutual repulsion [T51 US]- This 
basic feature of strongly-coupled Coulomb or Yukawa plas- 
mas is at the root of the quasi-localized charge approxi- 
mation, that has proved useful for the determination of 
dynamic quantities [20]: the charges inside the globule 
are trapped around local potential minima, and there- 
fore oscillate at a frequency close to the Einstein value 
w% = n a z 2 q 2 /(3me) [2I][22], where n Q = Z g / '(4tt R 3 g / '3) is 
the background density, and m is the counter-ion mass. 
The potential felt locally by a counter-ion then reads 



_ 1 z 2 q 2 n 2 _ 1 kTx 
1 ' ~ 2 3e 6 A 2 



(4) 



where x stands for the deviation from potential minimum, 
and A = [z 2 q 2 no/ '(efcT)] -1 / 2 can be thought of as a Debye 
length. The typical cage size is given by 6, as required from 
local electro-neutrality, and likewise, those cages located 
near the boundary of the globule (r = R g ) are centered at 
r — R g — 5 max - We expect <5 ma x and S to scale accordingly, 
and more precisely, we write <5 max = cxS, where a will be 
an important quantity for what follows. Due to the re- 
pulsion of neighboring counter-ions, we anticipate a < 1 
(the outer layer of confined ions is "pushed" towards the 
boundary r = R g , by a mechanism reminiscent of deple- 
tion in hard core systems). This is precisely the scenario 
at work in the ground state (i.e. at infinite T), where 
several approximate expressions have been proposed for a 
|15| . Consistent with these approximations and with nu- 
merical simulations that report 0.73 < a < 0.77 [TS], we 
will take a = 3/4 for our large coupling expansions. 

We are now in position to compute the uptake charge of 
the globule, from the number of counter-ions that arc able 
to escape their cage. At large V, only those cages located 
near the globule boundary can loose particles; there are 
(Rg/S) 2 such cages, so that 
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Since S 2 /X 2 = 3z 2 £b/S = 3T, and going from Z up to its 
rescaled form II, we obtain 



n oc r 3/2 exp(-a 2 r/2). 



(7) 



This shows, under strong coupling and at variance with 
mean-field, that the uptake charge actually decreases upon 
increasing the globule charge; furthermore, it is notewor- 
thy that our argument, valid at large T, also reproduces 
the mean-field small T behaviour with a power law of ex- 
ponent 3/2. 

We now seek a more microscopic information and at- 
tempt to predict the radial dependence of the counter-ion 
profile n(r) outside the globule. When a counter-ion ap- 
proaches the globule, it polarizes the trapped ions in their 
cages, and in turn feels the potential V(r) thereby created 
(the more obvious Z up /r contribution appears to be sub- 
dominant, see below). More specifically, the test particle 
located at r creates a field E = — zq(r — x)/(47re|r — x| 3 ) 
at point x where a counter-ion located inside the globule 
will be displaced from its equilibrium position, creating a 
dipole moment p = 3z 2 g 2 A 2 E/fcT [23]. The contribution 
of this dipole to the potential V felt by the test charge 
is p ■ (r — x)/(47re|r — x| 3 ), an expression that we have 
to integrate over all cages of counter-ions (one dipole for 
each cage); moreover, for consistency with our previous 
argument with outer cages centered at a radial position 

Ra 



V(r) 



<5max, and gathering expressions, we have 
3zVA 2 r 1 



-n 



d 3 x. 



(8) 
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To obtain the dominant contribution, in the vicinity of 
R g , we neglect the curvature of the globule, which gives 



-zqV(r) 



T 



R a 



with 



(9) 



The corresponding density profile follows from n(r) oc 
e -v/kT ^ -^rj^jQ^ yields the dominant behaviour logn oc 
(■b I {f — Rg + <5max)- This allows not only to compute the j2 
scaling parameter S but also to propose a scaling function 
for n(r). Indeed we obtain here S oc z 2 £ 2 3 /S 2 from the 
definition ©, i.e. S oc T 2 , and 
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n(r) 



oc F{Q = [(C + «)- 



(10) 



with £ = (r — R g )/S. As a consequence, if the scaling re- 
lation (flu)) holds, an important test for the consistency of 
our approach is to recover the same value of a as in expres- 
sion (JTj), close to 3/4. To summarize our analytical find- 
ings, we obtained, in addition to the explicit scaling form 
(JTUJ) , that while II oc S oc T 3 / 2 in the non-linear mean- 
field regime (meaning the limit of large bare charges Z g 
within Poisson-Boltzmann theory), strong coupling leads 
to S oc r 2 and an uptake charge (background plus counter- 
ions present within the globule) II oc r 3 / 2 cxp(— a 2 T/2). 
The slope (S) increases with T steeper than in mean-field, 
and ionic correlation effects lead at large V to a decrease 
of the charge II. 




Fig. 1: Counter- ion density profile n normalized by background 
density no as a function of distance from the globule center 
(normalized by the globule radius R g ). Here T = 13, R c /R g = 
2 and Z g — 200. The density peak location defines the distance 
5max, that is shown in the inset at various couplings, different 
cell radii (R c /R g = 2,4 and 8 and Z g from 50 to 3000. The 
dashed line shows the value 3/4, used for a throughout this 
work. 



To put these predictions to the test, we have performed 
Monte Carlo simulations, where the counter-ion interact 
through the exact Coulomb law and feel the background 
charge of the globule, the whole system being furthermore 
enclosed in a larger sphere of radius R c [2] . Following the 
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Fig. 2: Plots of rescaled counter-ion profiles in the vicinity 
of the globule, F" 1 log[n(r)/n(R g )], as a function of £ = (r — 
R g )/S. The different symbols correspond to different values of 
r, from 8 to 24. The thick dashed curve is for the function 
J-{Q/2 where T is defined in Eq. (|10|l . with a value a = 3/4. 
Inset: Integrated charge Qi Ilt vs radial distance, at weak and 
strong couplings. For T < 1, Qmt is always of the same sign 
as the background, while over-charging is observed at strong 
couplings. The dashed line shows the mean-field result. 



analytical treatment, we located the position of counter- 
ionic density peak close to the globule edge, which defines 
^max, see Fig. Q] The inset of this figure shows that the ra- 
tio a = S ma _ x /S does not depend on the coupling strength, 
and remains close to its ground state limit a ~ 3/4. Re- 
covering the proper value of a at large T can be viewed 
as an assessment of the validity of our simulations, and 
we can then proceed with the explicit check of the scaling 
form Eq. (|10p . It can be seen in Fig. [2] that the different 
curves exhibit good collapse at different V, and that the 
function F(C) captures the density decay in the external 
vicinity of the globule (£ > 0). While our argument above 
only provides the relation r _1 \og[n(r)/n(R g )] oc .F(C), a 
more refined analysis indicates that the prefactor is close 
to 1/2 [23] . so that we have plotted J r (C)/2 in Fig. H We 
have also computed the values of S and II in the simula- 
tions. They are shown in Fig. [31 which fully corroborates 
the analytical scaling behaviours. First, at small T, we 
have the mean- field behaviour S oc Ft oc T 3 / 2 , while at 
larger couplings, S oc T 2 and II becomes non monotonous. 
The detailed behaviour of II versus T provides a stringent 
test for our arguments: as shown in the inset of Fig. [3J the 
dependence of log(II/r 3 / 2 ) on T is linear at large T, with 
a (negative) slope compatible with the predicted value of 
ev 2 /2 = 9/32. We see that a unique value of a, inherited 
from ground state properties, accounts for the behaviour 
of the density profile together with more global quantities 
like the uptake charge. 

The previous considerations provide a detailed descrip- 
tion for the ionic density profile outside the globule. The 
behaviour for r < R g is more intricate, and has been an- 
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Fig. 3: Rescaled counter-ion density slope S defined in Eq. Q, 
and reduced uptake charge II = Z up z 4 £% /Rg , as a function of 
coupling parameter F. The dashed lines have slope 3/2 (to 
evidence the mean-field behaviour) and 2. Here S has been 
divided by 2 to have S ~ IT at low F. The inset shows II/r 3/2 
as a function of coupling, on a linear-log scale. The line has 
slope -a 2 / 2 - -0.28. 



alyzed in Monte Carlo. A quantity of interest is the total 
integrated charge Q- m t(r) inside a sphere of radius r. By 
definition, Q m t(R g ) = —Z up and Qint(-Rc) = due to 
electro-neutrality. Within mean-field, it is interesting to 
note that Qi n t is always of the same sign as the back- 
ground (here, negative). One can also note that the salt- 
free nature of our system imposes Q ln t{R g ) < 0, and more 
precisely, that Qi n t(r) < for all r > R g . Hence, a true 
overcharging cannot be observed in our salt-free system 
[ini[l2]- However, for T > 1, the inset of Fig. [2] shows 
that Qmt(f ) exhibits a range of distances where it is posi- 
tive: we refer to such a possibility as a local over-charging, 
and we note that it is somehow reminiscent of its counter- 
part occurring for impermeable colloids [TUHH]. Indeed, 
as for impermeable colloids, it is possible to prove that 
local overcharging is precluded within mean-field theories 
24,25], so that when present, it is a manifestation of ionic 
correlations, that become prevalent at T > 1. The elec- 
trophorctic consequences of this over-charging effect arc 
unclear, and left for future study. Incidentally, the close 
agreement between Monte Carlo and mean-field results at 
small T (see inset of Fig. [5]) can be seen as assessing the 
validity of the numerical methods employed. 

To summarize, we investigated the screening properties 
of a uniformly charged spherical globule, neutralized by 
point counter-ions. Invoking a quasi-localized charge ar- 
gument at strong Coulombic coupling T, we obtained an- 
alytically the counter-ion density profile n(r) outside the 
globule, together with two more global quantities: one, 
denoted S, follows from n(r) and is its characteristic in- 
verse decay length in the globule edge vicinity (r = R g ); 
the second quantity, n, stands for the reduced total charge 
inside the globule and quantifies the counter-ion uptake. 



At small couplings, S and n coincide and scale as T 3 / 2 . 
Strong ionic correlations, on the other hand, were shown to 
lead to a departure of both quantities: the slope 5" becomes 
steeper (algebraic increase in V 2 ) and the counter-ion up- 
take is much more efficient, leading to a charge n that de- 
creases upon increasing T (which can be obtained increas- 
ing the globule bare charge Z g ): the algebraic mean- field 
increase turns into an exponential decrease, see Eq. ([7]l. 
Our scaling predictions are free of adjustable parameters, 
and make use of known ground state properties [15] (that 
fix the parameter a = 5 max /S). These predictions were 
corroborated by Monte Carlo simulations, that provide 
the exact static properties of our system at arbitrary T, 
and that furthermore revealed an over-charging effect that 
is absent within mean-field scenario. We have treated the 
phenomenon of caging at a rather simple level, that turns 
out to be sufficient to capture the interesting violations of 
mean-field, driven by ionic correlations. In addition, while 
the mean-field theory applies schematically for T < 1 and 
the strong-coupling arguments cover the range V > 10, 
a quantitative understanding of the cross-over region at 
moderate couplings presumably requires intermediate ap- 
proaches in the spirit of Refs. [2~T] . 

Our work opens interesting venues for future studies. 
First, our predictions can be tested experimentally, in the 
spirit of the experiments reported in |26j . For the corre- 
sponding hydrophobic polyelectrolytes, we estimate that 
r ~ z 2 at room temperature, where z is the valency of the 
counter-ions. Consequently, with trivalent ions, one has 
r ~ 9, which is the regime where mean- field no longer 
applies, and our strong-coupling predictions take over, see 
Fig. |3] Second, the phase behaviour of an ensemble of 
such globules is unknown, together with the effects of an 
added electrolyte. Finally, the response to external per- 
turbations, both static (sedimentation) or dynamic (elec- 
trophoresis), should provide a relevant ground to investi- 
gate the signature of strong Coulombic correlations. 
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